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Delaunay Recovery of Cosmic 
Density and Velocity Probes 

W.E. Schaap Sz R. van de Weygaert^] 

Abstract 

Optimally resolved one-dimensional density and velocity profiles 
through cosmological N-body simulations are constructed by means 
of the Voronoi-Delaunay tessellation reconstruction technique. In a 
fully self-adaptive fashion a strikingly detailed view of the density 
features and the corresponding cosmic motions is recovered. 
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In essence, N-body simulations of cosmic structure formation are sup- 
posed to represent a discrete sampling of underlying continuous density 
and dynamical fields. The recovery of the corresponding continuous fields 
is a less than trivial exercise. They are often distorted by manipulated, 
user-dependent and therefore biased reconstruction schemes. This makes 
it in particular cumbersome to deal self-consistently with the character- 
istically multi-scale hierarchical nature of cosmological density fields. As 
significant is the failure to recover crucial structural aspects of the salient 
and frequently sharply defined anisotropic - filamentary and wall-like - 
patterns in the cosmic matter distribution. 

Recently, Schaap & van de Weygaert H have developed a fully self- 
adaptive and unbiased method to reconstruct density and related dynamic 
fields from a discrete and in general nonuniformly sampled set of point loca- 
tions. It is based on the stochastic geometric concept of Voronoi/Delaunay 
tessellations and forms an elaboration on the formalism first proposed by 
Bernardeau & van de Weygaert for the case of assessing the statistical 
properties of cosmic velocity fields. 

The application of the method to a large 256 3 GIF N-body simulation 
(LCDM, 141.3/i -1 Mpc, courtesy: S. White) g, § provides a beautiful illus- 
tration of its sizeable promise. The top panel of fig. 1 presents the particle 
distribution in a slice through this simulation. The corresponding density 
field determined through the Delaunay technique is shown in the adjacent 
panel. Notice how much better than the saturated particle plot this den- 
sity field manages to elucidate the wealthy and detailed structural features 
present in this cosmic volume, superbly rendering its high density contrasts. 
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Figure 1.1. Slice through a GIF N-body simulation (top), the Delaunay recovered 
density (center) and density and velocity profiles along the central line (bottom). 
(We are grateful to S. White for initiating these calculations.) 

While the image of the density field already provides evidence of its op- 
eration, it is through the objective assessment of the density profile along 
the central axis of the slice that the success of the method is fully man- 
ifested (solid line, lower panel). Evidently, the Delaunay technique yields 
a faithful representation of the density field over an impressive dynamic 
range, encompassing gently varying and extended low-density regions as 
well as the high density contrasts found in compact objects, be it either 
condensed clumps or the flattened dimension(s) of filaments and walls. 
Even more compelling is the correlation with the corresponding velocity 
field along the same line (dashed line). Largely superseding the poor ve- 
locity resolution in the conventionally shotnoise-dominated void regions it 
succeeds in reproducing the matter depleting super-Hubble like peculiar ve- 
locity flows (e.g. void at w 123/i _1 Mpc). Even more striking are the sharp 
velocity transitions encountered at the location of high density peaks, in- 
dicating the large induced infall motions along various directions towards 
these features (e.g. the peak at w 7/i _1 Mpc). 
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